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The possibility of identifying some of Galactic gamma-ray sources as clusters of primordial black holes is discussed. The 
known scenarios of supermassive black hole formation indicate the multiple formation of lower-mass black holes. Our 
analysis demonstrates that due to Hawking evaporation the cluster of black holes with masses about 10"'^^ g could be 
observed as a gamma-ray source. The total mass of typical cluster is ^ lOAf0. Detailed calculations have been performed 
on the basis of specific model of primordial black hole formation. 
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Introduction 



T-H Among numerous unsolved problems of astrophysics, 
there are two seemingly independent problems: first, the 
r-|j formation of supermassive black holes in galactic nuclei [U 
l_|H m 131 in m El IZl IE] and, second, the existence of unidentified 
gamma-ray sources observed by the Compton Gamma- Ray 
^ Observatory (CGRO) [llinilll]. The Energetic Gamma 
^-iRay Experiment Telescope (EGRET) operating in the en- 
Q ergy range from 30 MeV to 30 GeV [HI [12] is the main in- 
strument of this observatory. The existence of such sources 
has also been confirmed on the Fermi Gamma-ray Space 
, ^, Telescope FERMI equipped with the Large Area Tele- 
scope (LAT) on board, which covers the energy range from 
20 MeV to 300 GeV [Mj. 

The total view of the celestial sphere in the gamma-ray 
range was obtained for the first time on the CGRO (1991- 
2000) with the EGRET; this view revealed a surprisingly 
large number of unidentified point-like gamma-ray sources 
(170 of 271), according to the 3EG catalog of gamma-ray 
sources [9]. However, owing to a certain improvement of 
the diffusion radiation model used to analyze the observa- 
tional data, the 3EG data were revised and the new EGR 
catalog of the sources was composed (TUj, where the num- 
ber of unidentified sources is 87 of 188 and the distribution 
of the sources over the celestial sphere is more isotropic. 

The Fermi Gamma-ray Space Telescope FERMI was 
successfully launched into the geostationary orbit in sum- 
mer 2008; the results obtained on this telescope were pre- 
sented in February 2010 in the first FERMI IFGL catalog 
[T3] . In contrast to the 3EG and EGR catalogs, the IFGL 
catalog contains 630 unidentified sources of 1451. Their 
distribution on the celestial sphere is much more isotropic 
that those in the 3EG and EGR catalogs. 
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The second problem of modern astrophysics is the ori- 
gin of supermassive black holes (SMBH). Various forma- 
tion scenarios of massive primordial black holes were dis- 
cussed in [2 [21 HH US] . Below we briefiy enumerate list the 
proposed mechanisms of primordial black holes (PBHs) 
formation (see also the detailed review pT)). 

PBHs were first suggested as a result of adiabatic fluc- 
tuations at the radiation-dominated stage of evolution of 
the Universe [H [HHnHII] ■ Masses of such black hole do 
not exceed solar mass. 

Variety of models with flat or bumped spectra of isocur- 
vature fluctuations were proposed [2 ^ (23 1 (24 1 [25 l (26 1 (27 1 
[2S1 Uni 130] • A very promising mechanism of SMBH forma- 
tion from the large amplitude isothermal fluctuations in 
baryonic charge density was proposed by Dolgov and Silk 
|31j . In this mechanism the bumped spectrum of isother- 
mal fluctuations is a byproduct of vacuum bubble collapses 
during phase transition at the inflationary stage. Masses 
of black holes could be larger the solar mass in orders of 
magnitude in this case. 

Spectrum of fluctuations depends on the form of the 
inflaton potential, see e. g. [32l[33]. As was shown in [34] . 
it opens possibility to produce high density fluctuations 
that collapse afterwards into SMBHs. In this case, masses 
of black holes could be much larger than the solar mass. 

Black holes with masses ~ M0 possibly formed at 
quark-hadron phase transition at the cosmological time 
10~^ s, [35l [35. Today such black holes would be a com- 
ponent of dark matter . 

Suppose that some massive non-relativistic particles 
dominates at an early stage of the Universe evolution. In 
that case the pressure is negligible and could not prevent 
gravitational collapse of high density regions [21 . 

Usually, false vacuum decay is accompanied by forma- 
tion of spherical walls with a true vacuum inside [371I551I551 
[ini SD [11 [H [M [IS] . The walls quickly expand and col- 
lide what could lead to collapse of islands of false vacuum 
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into black hole. Most massive black holes of order IMq are 
produced during the period of quark-hadron phase transi- 
tion. 

The mechanism of black holes formation from the closed 
domain walls was proposed in |16l I46j and developed in 
|47) . These domain walls could be originated due to evo- 
lution of a scalar field during inflation. An initial non- 
equilibrium distribution of scalar field imposed by the back- 
ground de-Sitter fluctuations gives rise to the spectrum of 
PBHs, which covers a wide range of mass — from small 
masses up to supermassive ones. The PBHs of smaller 
masses are concentrated around the most massive ones 
forming a fractal-like cluster. It was revealed that this 
mechanism is a rather common for many inflationary mod- 
els and it is worth to discuss it. 

The cosmological models of the formation of super- 
massive primordial black holes can hardly be free of the 
multiple formation of less massive black holes (see, e.g., 
[ini SH] ) ■ Moreover, to suppress the predicted abundance 
of black holes of smaller mass, special assumptions are re- 
quired [49^. We devote special attention to the model of 
PBH formation in the next sections. 

The detection of black holes with masses below IO^Mq 
is difficult in view of the weak accretion of surrounding 
matter. Another method for detecting black holes by Hawk- 
ing radiation [50] is efficient only in the vicinity of the 
Earth. Indeed, even in case of homogeneous distribution 
of PBHs over Galaxy, the rate of single PBH observation 
would be ~ 10^^ years" ^. Meantime clusters containing 
a large number of black holes with masses below or about 
lO^^g are visible from larger distances. 

In this work, we demonstrate that clusters of primor- 
dial black holes can be detected as point-like gamma sources 
due to effect of Hawking evaporation. According to our es- 
timations such a cluster contains black holes with different 
masses starting from tens of solar mass and smaller. Num- 
ber of small black holes appears to be large enough so that 
their total radiation could be seen from the Earth. As a 
result, the origin of unidentified gamma-ray sources can 
be explained and the way for detecting massive primordial 
black holes is indicated. 

2. Clusters of black holes from closed domain walls 

In this Section we briefiy discuss specific model of mas- 
sive black hole formation. The details can be found in pa- 
pers [ini 27] where the mechanism of massive primordial 
black holes production was elaborated. 

Some infiationary models suppose a creation of our 
Universe either near a maximum of potential of inflaton 
field or near its saddle point (s) to realize a desired slow 
rolling providing a sufficient number of e-folds (see details, 
e. g. in [5TJIS2]). As it will be shown below these models 
include the possibility of the formation of macroscopically 
large closed walls from a scalar field. After the end of in- 
flation these closed walls collapse to BHs if these walls are 



large and heavy enough [Tni US]. This mechanism is re- 
alized in well known models like the Hybrid Inflation [S3] 
and the Natural Inflation [M]. A scalar fleld could be the 
inflaton itself or some additional field. 

Consider general mechanism of closed wall formation 
based on quantum fiuctuations near unstable point (s) like 
a saddle point or a maximum of potential of scalar field. 
An evolving scalar field may be split into classical part, 
governed by the classical equation of motion, and quan- 
tum fluctuations |55j . To facilitate the analysis, let us 
approximate the potential near its maximum as 



V = V„ 



m 



(1) 



where without the loss of generality the maximum is as- 
sumed at = . Then, the probability density to flnd a 
certain fleld value </> has form [SB] (adapted to the consid- 
ered case): 



dP(0,T;0i„,O) =d0 



X exp 



(2) 



Here a ~ ^/cr^, ^ = m'^/3H and a = ff'^/^/27r, where the 
Hubble parameter H ~ ,/87rVo/(3Mpi). 

Let us choose a positive value for the initial field, 0in > 
0. Then an average field value will increase with time, 
ultimately reaching the minimum of the potential at some 
value > 0. This means that a greater part of space will 
be finally filled with the field value = Meanwhile, 
the field in some (small) space domain could jump with 
the probability ^ over the maximum due to the quantum 
fluctuations. In the following, an average value of the field 
representing this fluctuation tends to an another minimum 
of the potential, 0_ < 0. As the result, space at the final 
stage will be filled by vacuum while some space domain 
is characterized by the field value = i)f)_ < 0. If one starts 
to move from inside of the domain to the outside, the path 
would start from a space point with (/>_ and finish at a 
space point with (j)^. Hence, the path must contain the 
point with the maximum value of potential. It means that 
a wall is formed inevitably between such space domains 
and the "outer" space with (p = (/)+ [56j [57] . 

The "dangerous" values of fluctuations are those with 
4> < 0. Such space domains will be surrounded by closed 
walls and if their number is sufficiently large it would 
strongly infiuence the dynamics of the early Universe. If 
a fraction of space surrounded by the walls is not very 
large, the resulting massive BHs, which are formed from 
the walls, could explain the early formation of quasars [58] . 

3. Specific features of black hole clusters 

We will be founded on the described above scenario 
of primordial black hole formation, see also [S] H] dni 
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B71 [55] , This scenario predicts the existence of super- 
massive black holes in galactic nuclei and the existence 
of intermediate-mass black holes in galactic halos at large 
distances from their centers. One of the feature of this sce- 
nario is the formation of black holes with a characteristic 
cluster structure. 

The scenario of the formation of massive primordial 
black holes cannot predict their present mass distribution, 
since it strongly depends on initial conditions in the pe- 
riod of inflation and the parameters of initial Lagrangian. 
Moreover, the initial distribution of primordial black holes 
is distorted during the subsequent evolution of a galaxy 
when clusters of black holes are merged with each other 
and with the supermassive black hole in the galactic center. 
The remaining primordial black holes form a population 
of black holes in the galactic halo. Thus, for our aims, it is 
sufhcient to determine the initial conditions under which 
the number of formed primordial black holes is certainly 
larger than the number of the unidentified gamma sources. 

The procedure for determining the mass spectrum of 
primordial black holes was considered, e.g., in [2l|6], where 
the potential 



-.9) (3) 



was chosen for certainty. It was shown that at the rea- 
sonable magnitudes of parameters A, / and A (/ = 10.0 
and A = 1.66 in the units of Hubble parameter at the in- 
flationary stage and almost arbitrary A), one obtains an 
appropriate structure of clusters. PBHs within a cluster 
have mass distribution approximated by the power law (see 
Fig. [T]) with typical total mass ^ IOMq. 

/. v. yc 1 1 ■ / n/i \ 

= /i„(Afi„) 



dN 
dM 



Al 



(4) 

Black holes with masses Mi„ ~ Af„ ^ 10^^ g are of partic- 
ular interest, because they exist at the final, most intense 
evaporation stage now and produce gamma-ray fluxes. More- 
over, observational constraints on them are particularly 
strong [15]. 

As was mentioned above, the resulting mass spectrum 
is strongly distorted due to merging of black holes, ac- 
cretion, and evaporation. The last process is important 
for low-mass black holes with masses M < Af* ^ lO^^g. 
Note that non-vanishing mass distribution in this region 
is supported by Hawking evaporation of bigger black holes 
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Figure 1: Mass distribution of black holes in a cluster, x = M/Mt. 



the energy of the particles 
0- = 277rA/Vm|i. 



(6) 



As a result, the black hole mass loss rate is given by the 
simple expression 



dM 
~dt 



dE-y 
dt 



const 



(7) 



where the coefficient k presents the contributions from all 
of the particles divided by the contribution from photons. 
Neglecting the dependence of k on Af and using Eq. 0, 

we 

obtain the known relation between the initial and present 
values of the mass of black holes: 



Mi„ = (Af3 + A/; 



,1/3 



(8) 



In this case, the mass distribution of black holes in the 
cluster given by Eq. ^ is transformed as (see Fig. [T]) 



fiM) 



dAfir 

dM 



fin{Min{M)) = 



4.4 X 10 



17 



{M/M^ 



Af* 



1 + {M/M.,f 



4/3 ■ 



(9) 



The predicted total number of black holes in the cluster 
and its mass in the present Universe are N = J f{M)dM ~ 
4 X 10^^ and Aftot ~ / Mf{M)dM ~ 9.5Mq, respectively, 
within the same cluster. Let us estimate the present evaporation- The initial conditions of the formation of black holes in 



induced distortion of the mass spectrum of primordial black 
holes. 

The black-hole evaporation rate is characterized by the 
temperature |50] 



T ■ 



1 m 



PI 



21^ MeV, 



(5) 



Stt M M 

where mpi = 1.2 x 10^^ GeV is the Planck mass. 

In the approximation of high radiation energies, the 
effective evaporation area of black holes is independent of 



the simulation are chosen in such way that the number of 
the formed clusters of black holes is certainly larger than 
the number of gamma-ray sources. Under the chosen con- 
ditions, the galactic halo in the early Universe contained 
about Nc\ ^ 1400 clusters with the typical sizes -R < 1 pc. 

The lifetime of such a cluster due to black hole escape 
can be estimated as [3 [SH] 



t ~ 40t,.ci, 



(10) 



3 



where t^ci is the relaxation time [51 [SS] 



1 

4tt G^m'^nlogOAN' 



(11) 



4. Photon signal from clusters of black holes 

Let us estimate the luminosity of the typical cluster 
of PBHs. The intensity of the direct emission of photons 
from one black hole with the temperature T is given by 
the expression [15] 



~ 2^exp{E/T) 



dE dt 



(12) 



and strongly depends on the mass of the black hole. The 
convolution of spectrum ( 12 ) and mass distribution ^ 
gives the rate of photon emission from the cluster of black 
holes (at E > 100 MeV) 



TV « 6.6 X 10^^ s"\ 



(13) 



The total luminosity of the cluster of black holes is calcu- 
lated as 



dN 

L(M) dM 

^ 'dM 



9.7 X 10'^'^ GeV/s, (14) 



where the luminosity of single black hole is taken into ac- 
count in the form 50, 



= 2.2 X 10^ 



A/, 
If 



GeV/s. (15) 



Note that the mean energy of emitted photons is i5 ~ 
Lci/N 15 MeV and, hence, consideration of clusters of 
black holes as gamma-ray sources is justified. 

Let us determine the number of clusters of PBHs con- 
sidered as gamma-ray sources that can be observed on the 
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V ~ GmN/R ~ 3 X 10'^ cm/s is the velocity of black O 
holes in the cluster and m, N, and n ^ N/W are mass, 
number, and density of typical black holes in the cluster ^ 
respectively. At the end of this interval the gravitational 
collapse takes place. Detailed discussion applied to this is- 
sue is given in [5] . The lifetime appears to be many orders 
of magnitude larger than the age of the Universe; there- 
fore, it can be thought that the total number of clusters 
survive to the present epoch. Despite a considerably low 
mass of the cluster, its lifetime is long because of the high 
density of black holes in the cluster. However, no exter- 
nal destructive influences were taken into account so the 
fraction of survived clusters is quite uncertain parameter. 

The strongest limit on the density of PBH relates to 
mass interval around ^ 10^^ g [15] . The density in the 
present Universe ripBH ^ 3 x , obtained in our cal- 

culation, does not contradict this limit. 
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Figure 2: The expected photon spectra from PBH clusters (the 
"banded" region). Below the region the flux is under LAT integral 
sensitivity (r > -Rmax), and the expected number of clusters is less 
than unity above the region. The differential sensitivity of detectors 
INTEGRAL and OSSE are also shown. 



LAT. The photon flux arriving at the Earth from such a 
cluster is 



E =^ 



(16) 



which must exceed, in order to be detected by Fermi LAT, 
the threshold value [T^ 



i^"™ = 3 X IQ-^ cm-^s-i. 



(17) 



Then, the maximum distance i?inax at which the radiation 
source can be resolved is 



1 N 

47r Eru„ 



1/2 



4.3 X 10^ pc 



and therefore the number of such sources is 

3 



Tici X -7ri?„ 



R 



gal 



cl 



33, 



(18) 



(19) 



where rid is the number density of clusters, which are as- 
sumed to be located within the characteristic size of the 
galaxy i?gai 15 kpc. 

The gamma-spectrum seems to be a promising tool for 
PBH cluster identification. For energy of interest {E > 
100 MeV) Eq. (|9| and Eq. ^ give 

In the fig.[2]we present the fluxes from PBH clusters. Note 
that the maximum at i? ~ 10 MeV in the predicted spec- 
trum corresponds to the peak of mass spectrum fig. [T] at 
Mq ~ IQi^ g. 

As seen from fig.[2]X-ray telescopes are unable to detect 
PBH clusters because their sensitivity is much less than 
the threshold value [BUI [H] • 
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Figure 3: Unidentified Gamma ray sources seen by Fermi LAT (blue 
cross) on celestial sphere are shown. The red filled markers are 
sources with spectrum indices 3 within la error, unfilled markers 
are sources with spectrum indices 3 within 3cr error. 



The Fermi LAT has detected 15 sources with spectral 
index — 3 at Icr and 93 sources at 3cr error [13 . Their distri- 
bution on the celestial sphere is isotropic and their amount 



is in good agreement with the prediction of Eq. ( 19 ) (see 

Each PBH cluster should be also a point-like source of 
high-energy neutrinos. We estimated the expected neu- 
trino flux in similar manner and compared it with sen- 
sitivity of neutrino observatories. It was obtained that 
neutrino observatory AMANDA [62l [63], which is sensi- 
tive to the flux of neutrinos with energy > 1.9 GeV ex- 
ceeding ^ lO^^^cm^^s^^, can potentially detect neutrinos 
from the cluster at several kpc. At these distances we can 
expect about one cluster. It means that neutrino experi- 
ments (AMANDA, IceCube) are promising for indication 
of PBH clusters. 

5. Conclusions 

In this work, we attempted to unify two astrophysical 
problems. First, we propose the explanation of the ori- 
gin of unidentified gamma-ray sources. Second, new tool 
of low-mass black holes searching is discussed. Multiple 
production of the latter is predicted by some models of 
the formation of massive primordial black holes in galactic 
centers. It has been shown that if the spatial distribution 
of PBHs has a cluster structure, an individual cluster of 
PBHs is detected as a point-like gamma-source. Total lu- 
minosity of small black holes in a distant cluster is large 
enough to be detected on the Earth. Note that the model 
of the formation of massive primordial black holes pro- 
posed in [ini SS] predicts just the cluster structure. The 
mass spectrum of PBHs, their abundance, and their radia- 
tion energy spectrum calculated using this model indicate 



the real possibility of the correlation between point-like 
gamma-ray sources and PBHs. 
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